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Abstract 

We calculate the electromagnetic corrections needed to obtain isospin invariant 
hadronic pion-pion s-wave scattering lengths , from the elements Occ^a^c of 
the s-wave scattering matrix for the (tt^vt^, tt'^tt'') system at the tt+tt" threshold. 
The latter can be extracted from experiments on the vr+vr" atom (pionium). Our 
calculation uses energy independent hadronic pion-pion potentials , V'^ that satis- 
factorily reproduce the low-energy phase shifts given by two-loop chiral perturbation 
theory, with the hadronic mass of the pion taken first as the charged pion mass and 
then as the neutral pion mass. We also take into account an important relativistic 
effect whose inclusion influences the corrections considerably. 
PACS: 13.75.Gx,25.80.Dj 
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1 Introduction 



The extraction of the s-wave pion-pion scattering lengths , J = 0, 2, from the 
results of the DIRAC experiment currently in progress at CERN [1] requires 
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a knowledge of the electromagnetic corrections oqc — ^(a^ — a°) and — 
(|a*^ + |a^), where Occ and aoc are elements of the well known matrix K 
of scattering theory in the s-wave for the two-channel (tt+tt^, 7r°7r°) system 
at the TT+TT" threshold. (The subscript c refers to the n^n^ channel, to 
the 7r°7r° channel.) Previous attempts to calculate these corrections using a 
potential model have been reported in Refs.[2,3]. In this paper we improve on 
the calculation of Ref. [3] by using new hadronic potentials obtained by fitting 
the phase shifts given by two-loop chiral perturbation theory (ChPT) and by 
including a relativistic modification which turns out to have a significant effect 
even in the threshold situation we are considering. 

This relativistic modification emerged as an important issue in the calculation 
of electromagnetic corrections for the analysis of low energy 7r~p scattering 
data [4] . The treatment of the two-channel (7r~p, 7r°n) system closely parallels 
that of the (tt'^tt", 7r°7r°) system and the coupled relativised Schrodinger equa- 
tions (RSEs) for the two systems are formally identical. The two-channel RSEs 
that we now use to model the s-wave of the physical (Tr"*"?!", 7r°7r°) system are 
therefore given by Eq.(21) of Ref. [4]: 

l2y^ + Q' - 2mfV(r) U(r) = 0, (1) 



where 



In Eq.(2), /ic and /iq are the physical masses of tt^ and tt^ respectively, while 
Qc and qo are the cm. momenta of the two channels. Further, 

_ W - 2^1 _ W - 2^1 



where W is the total energy in the cm. frame. Note that fc and /o are 1 at 
the respective thresholds W — 2//c and W — 2//o and increase as W increases. 

Eq.(l) now replaces Eq.(8) of Ref. [3]. The inclusion of the factor f is of crucial 
importance for the calculation of the electromagnetic corrections for both the 
(tt+tt", 7r°7r°) and (7r~p, 7r°n) systems. The first part of Section 3 of Ref. [4] 
gives the arguments for its inclusion in the RSEs for the latter system. The 
potential matrix V has the form 

Y ^ yem _^ yh_ 
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The electromagnetic potential matrix has only one nonzero entry (V^™')tc = 
V^"^. In the notation of Ref. [4], K"^™ contains only V^"^ (point charge Coulomb) 
and V^^* (which takes account of the extended charge distributions); the effect 
of y and V'"^ (vacuum polarisation) on the positions and widths of the levels 
of pionium is treated separately. 

The hadronic potential matrix V'* is assumed to be isospin invariant, in accor- 
dance with the arguments of Gasser and Leutwyler [5,6], which show that there 
is practically no isospin breaking of the purely hadronic pion-pion scattering 
amphtudes. Thus 



The determination of the hadronic potentials will be considered in Section 2. 
Once they are chosen, the coupled equations (1) are integrated numerically to 
give K as a function of W. Extrapolation to W — 2ijLc gives the matrix 



whose elements ttcc, doc can be obtained from experiments on pionium. The 
determination of these quantities from pionium data is reviewed at the be- 
ginning of Section 3, which then sets out the results of our calculation of the 
electromagnetic corrections needed to obtain the hadronic pion-pion scatter- 
ing lengths a°,a^ from Occ and Oqc- In Section 4 we compare our results with 
those of ChPT. 



2 The hadronic pion-pion potentials 

We now need to model the hadronic situation and determine the potentials V^, 
V'^. Chiral perturbation theory (ChPT) has proved a fruitful effective theory 

of low energy pion-pion scattering that incorporates the essential constraint 
of chiral invariance. From the point of view of ChPT, the hadronic mass fi 
of the pion (assumed the same for tt^ and 7r°) is an adjustable parameter 
and ChPT is able to generate hadronic phase shifts 5°, 5"^ that depend on the 
value of II that is chosen. This procedure is complicated by the need to fix 
certain low energy constants using experimental data obtained with physical 
pions. Nevertheless, Colangelo [7] has provided us with hadronic s-wave phase 
shifts d^{iJ,; W), I = 0,2, obtained from two- loop ChPT as functions of W for 
a region above the threshold W — 2/1, for /j, = and for = /iq. 





a = K(2/.,), 



(6) 
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Table 1 

The hadronic pion-pion scattering lengths for total isospin I = 0,2 and hadronic 
mass = //c, A*o given by two-loop ChPT (Ref.[7]) and by the potentials F° and 
of Fig.l (all results in fm). 





two-loop ChPT 


our potentials 


a°(Mc) 


0.3055 


0.3004 


a°(Mo) 


0.2930 


0.2892 




-0.0632 


-0.0632 


a^(Mo) 


-0.0614 


-0.0613 



What we have been able to do is to find potentials V^, V'^ that are independent 
of both the energy W and the hadronic mass /x, which reproduce quite well 
the two- loop ChPT phase shifts of Ref . [7] for /i — /ic and /i — /iq via the RSEs 

{^2 + - l^cfcV\r)^ u{r) = 0, (7) 
^2+^0- y^^hV\r) \ u{r) = 0, (8) 



with / = 0, 2. The important difference from Eq.(6) of Ref. [3] (which applies 
to 11 = lie only) is the presence of the factors fc and /o in Eqs.(7) and (8). 
Compared with the work of Ref. [3], we were also able to use our experience 
of constructing potentials for calculating electromagnetic corrections in low 
energy pion-nucleon scattering to obtain physically more realistic potentials 
than the double square wells used earlier. 

We found that the best fit to the phase shifts was obtained when the range 
parameter in these potentials was 1.5 fm. Each potential contains three fur- 
ther parameters that were varied to obtain the best fit. Full details of the 
parameterisation are given in Ref. [8]. We found that it is possible to construct 
for / = 0, 2 potentials V\r) that give a quite good fit to the two-loop ChPT 
phase shifts 5^{ii] W) for ii — jic and ji — jiQ, using Eqs.(7) and (8) to gener- 
ate phase shifts by integrating the regular solution outwards from the origin. 
These potentials are plotted in Fig.l. It is impossible to obtain a satisfactory 
fit to the phase shifts for both values of /x when the factors fc and /o are 
omitted from Eqs.(7) and (8); this is further significant evidence for the need 
to include them. What we have been able to do is to capture some of the 
dynamical behaviour of ChPT by means of a simple potential model involving 
potentials independent of the energy and the hadronic mass. In Figs. 2 and 3 
we show how well the potentials reproduce the ChPT phase shifts for fi = fi^ 
and fi = fj,o respectively. Table 1 gives a comparison of the hadronic scattering 
lengths given by two-loop ChPT and by our potentials. 
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Fig. 1. The potentials V'^{r) and V'^{r) that give the best fit to the two- loop ChPT 
phase shifts of Ref.[7] for hadronic masses /Xc and f^o. 

3 Prom pionium data to hadronic pion-pion scattering lengths 



The first step is the extraction of the quantities Ucc, ooc from the experimental 
values of the lifetime r of the Is level of pionium and of the difference AW^"''^ 
between the energies of the 2s and 2p levels. This is considered in Section 2 
of Ref . [3] , with the results 



aoc(fm) = -0.3794 [r(eV)]^/2^ 



(9) 



acc(fm) = -0.4167 AW''"^(eV) 



(10) 
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Fig. 2. The phase shifts 5° {fie, W) and S'^inc, W) (in de grees) for hadronic mass //c 
given by two-loop ChPT (crosses) and by the potentials y° and of Fig.l (solid 
curves). 

In terms of the lifetime r, Eq.(9) becomes 

aoc(fm) = -0.3078 [r(fs)]-^/^ (11) 



The slight differences between the numerical constants of Eqs.(9) and (10) 
and those in Ref.[3] come from a new value of Occ and the direct inclusion 
of the effect of vacuum polarisation in Eq.(9). These results are based on the 
work of Ref. [9], which uses general principles of scattering theory, in particular 
the analytic continuation of K below the tt+tt^ threshold. The basic formal- 
ism goes back to Hamilton, 0verbo and Tromborg [10] and was worked out 
systematically for multichannel situations by Rasche and Woolcock [11]. It is 
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Fig. 3. The phase shifts (5°(^o; W) and 5'^{no; W) (in de grees) for hadronic mass //q 
given by two-loop ChPT (crosses) and by the potentials and of Fig.l (solid 
curves). 

independent of the potential model, which we use only to calculate electro- 
magnetic corrections. 

It is important to observe that the expressions for the hadronic shifts and 
widths of the levels of hadronic atoms obtained in Ref. [9] are given in terms of 
the elements of a as defined in Eq.(6). As shown in Ref. [9] there is no signif- 
icant error involved in evaluating K at the threshold of the charged channel 
instead of at the pole position corresponding to the bound state. Thus the 
extraction of Occ and aoc from experimental data is completely separated from 
the question of the electromagnetic corrections that relate these two num- 
bers to the hadronic s-wave scattering lengths , I = 0,2. This separation is 
valuable from a conceptual point of view. 
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In order to compare our results with field theoretical treatments of the width 
r, it is necessary to give our definitions of two-channel scattering amplitudes 
and to establish the connection with the standard amplitudes used in work on 
ChPT. We confine ourselves again to the s-wave. Starting from the unitary 
matrix S, we define matrices T, K and A as follows: 



T = -z(S-l2)(S + l2)-\ 

Ql/2rj.-lQl 



Oil \^ 1 



+ 



-Ph{rj^) 




[12) 



, (13) 



(14) 



The matrix Q is given in Eq.(2), while /3 = a/ic, V- — ~P/'^Qc, 



Eq.(13) is the explicit definition of K, in which the specific low energy be- 
haviour of T induced by the Coulomb interaction is taken into account. The 
complex quantity A^c plays the key role in the formalism of Ref . [9] . 



If we denote by A^ the s-wave part of the matrix of scattering amplitudes 
used in work on ChPT, the connection with the matrix A of Eq.(14) is 



A = (327rW)~^ A^ 

1/ I 1 



(15) 



Note that A has the dimension of a length, while A^ is dimensionless. The 
most satisfactory field theoretical treatment of pionium is given in Gall et a/. [12], 
which uses QCD (including photons) and works with the low energy expansion 
provided by ChPT. Ref. [12] also gives references to and discusses other field 
theoretic methods of treating hadronic atoms. In Ref. [13], Gasser, Lyubovit- 
skij and Rusetsky apply the formalism of Ref. [12] to the numerical calculation 
of r. A more detailed account of their work is given in Ref. [14]. The expression 
for r given in Refs. [12-14] involves the real part of an amplitude ^4^"°*^, from 
which the specific low energy behaviour due to the Coulomb interaction has 
been removed. Using Eq.(15), we have 

Aoc(2/Xe) = y2(647r/Xe)"'A+;°°- (16) 



If we compare the expression for T given in Refs. [12-14] with our result in 
Eq.(9) of Ref. [3], on which Eq.(9) above is based, we find that a term pro- 
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portional to In a and the purely hadronic scattering lengths appear in the 
quantity K given in Eq.(4) of Ref.[12], while our result has no term involving 
In a and contains the scattering parameters that include the effect of electro- 
magnetic interaction. As Ivanov et al. [15] have pointed out, these differences 
come from the same source: the treatment of Coulomb photons. In our method 
the Coulomb potential is included in the RSEs (1), which are solved exactly. 
The nonanalytic behaviour in a is implicitly included in our electromagnetic 
corrections. The treatments based on scattering theory [9] on the one hand 
and on field theory [12-14] on the other therefore lead in the end numerically 
to the same results. 

To understand the electromagnetic corrections we must distinguish between 
the physical situation, in which the electromagnetic interaction is present and 
the pions have their observed masses, and the hadronic situation, in which 
the electromagnetic interaction is switched off. The pioneering work of Gasser 
and Leutwyler [5,6], already used in Section 1, shows that, at first order in 
the low energy expansion of ChPT, in the hadronic situation the only isospin 
breaking effect due to the difference between the quark masses rriu and rud is 
a very small shift in the 7r° mass. It is therefore a very good approximation to 
take the hadronic situation as one in which all the pions have the same mass 
II and the s-wave amplitudes are isospin invariant. 

As noted in Section 2, for ChPT at the hadronic level, the value of is a 
parameter that is not determined by the theory. (This does not mean that its 
value can be chosen completely arbitrarily.) Thus ChPT contains dynamical 
information about s-wave hadronic tttt scattering for a whole range of hadronic 
masses of the pions. The potentials constructed in Section 2 capture a lot 
of this dynamical information for a range of masses around Hq. They were 
constructed in order to be able to incorporate the effects of the electromagnetic 
interaction (the Coulomb potential for extended charges and the difference 
between /xq and /Xc) into a potential model of the physical situation that enables 
the electromagnetic corrections to be calculated. That model is just the two- 
channel RSEs(l); it is simple and hopefully will give a reasonably reliable 
estimate of the corrections. 

We now define the hadronic quantities a'^cil^), o,Qcil^) by 

a'M = + la'ili) , aHii) = ^{a\ii) - a\ii)). (17) 

As stated in Section 3 of Ref. [3], numerical experience shows that the quanti- 
ties least sensitive to small variations of the hadronic potentials are 

Aaccifi) = a^c(/i) - ttcc , Aoc(/^) = aQ^{n)/aoc - 1. (18) 
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The term 'electromagnetic corrections' will in future mean Aacc(/i) and Aoc(/u). 
Note that the hadronic quantities a\fi),I = 0,2, dccif^) and aQ^(/i) depend 
exphcitly on /x. They are given by an equation like Eq.(7) or Eq.(8), with 
replacing or /xq and a momentum variable q{n; W) corresponding to the 
mass /I. However, the physical scattering parameters Ucc and qqc are given by 
the RSEs(l) and do not depend on fi. The a^{fi) given by the potentials of 



Section 2 are listed in the second column of Table 1 for ji — jic and ji — jiQ. 
They lead to the results 

a!l^{Hc) = 0.1791 fm , a^c(/^c) = -0.1714 fm, (19) 

a^(//o) = 0.1721 fm , = -0.1651 fm. (20) 



The procedure for determining the elements of K at the threshold of the 
charged channel was described briefly at the end of Section 1. The coupled 

RSEs(l) were integrated numerically to obtain K as a function of W , as 
described in Section 3 of Ref. [16]. The matrix K was then extrapolated to 



W = 2jjic to obtain a. The results are: 

Occ = 0.1773fm, aoc = -0.1726fm, aoo = 0.0576fm. (21) 

Prom Eqs.(19)-(21) the electromagnetic corrections of Eq.(18) are then 

Aacc(/^c) = O.OOlSfm, Aoc(/ic) = -0.0068, (22) 
Aacc(Ato) = -0.0052fm, Aoc(Ato) = -0.0436. (23) 



Eqs. (22) and (23) are the final results of our calculations of the electromagnetic 
corrections in the potential model. They supersede the results in Refs.[2,3], 
where much more primitive hadronic potentials were used and the relativistic 
factors /c, /o were not taken into account. 

We have given the corrections for two possible values of the hadronic mass 
of the pion. However it was also shown by Gasser and Leutwyler [5,6] that 
the mass of the pion in the hadronic situation is very close to //q- It follows 
that the hadronic starting point, with respect to which the electromagnetic 
corrections need to be calculated, has a pion mass /io- The results in Eq.(23) 
are therefore the true electromagnetic corrections. These need always to be 
calculated with //q as the hadronic mass of the pion. 
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4 Discussion 



We now compare our results with those obtained from the work of Knecht 
and Urech [17], who use the low energy expansion of ChPT and include the 
effect of the electromagnetic interaction. At lowest order in ChPT the hadronic 
scattering lengths are 



givmg 



""cM = ^ ' «oc(/^o) = > (25) 



where F is the pion decay constant in the chiral limit (~ 88 MeV). Using 
Eq.(15) above and Eqs.(2.14)-(2.16) of Ref.[17], we find that 

_ (2^g-/ig) 1 _ V2(4/xg-/ig) 1 



for the physical situation, with the pions having their observed masses. Com- 
bining these results, the lowest order calculation of Ref.[17] gives the correc- 
tions 

Aacc(Ato) = -0.0139fm, Aoc(Ato) = -0.0533. (27) 



At the one-loop level, only the amphtude for n'^n — > n^n^ is calculated in 
Ref.[17]. Using Eq.(15) again, the results in Eqs.(5.2) and (5.6) of Ref.[17] 
yield the quantities 

a^^iHo) = -0.1570 fm , oqc = -0.1675 fm, (28) 



from which we deduce that 

Aoc(/^o) = -0.0625, (29) 



to be compared with the leading order result in Eq.(27). These numbers have 
been calculated using the values of the low energy constants = 1 — 4, and 
/Cf°, JCt^ given in Ref.[17]. 
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To compare the results in Eqs.(23) (present work) and (29) (calculated from 
Ref. [17]) we separate the effect of the pion mass difference from that of the 
Coulomb interaction. The correction due to the mass difference only is 



A^fiHo) = -0.0443 (present work), A'SfifXo) = -0.0561 (Ref.[17]) . (30) 



The uncertainty in the second number, based on the quoted uncertainties in 
= 1 — 4, is 0.0015. Moreover, the change from the lowest order result to 
the one-loop result is small. The number from our present work is subject 
to an uncertainty that comes from the possible difference between the purely 
hadronic potential that we have used and the effective hadronic potential in 
the presence of the electromagnetic interaction. This difference can be thought 
of in a field theory picture as the result of electromagnetic mass differences 
in the internal lines of Feynman diagrams. Assuming that these effects are 
correctly accounted for in Ref. [17], it is likely that the ChPT result gives a 
more reliable estimate of the mass difference correction than does our potential 
model. 

On the other hand, due to the large uncertainties in ICf^, A^2 °) the Coulomb 
correction in the ChPT calculation has a large uncertainty; it is —0.0064(82). 
The potential model gives a small correction, +0.0007, with no significant 
uncertainty. This result lies within the ChPT error band and removes the 
uncertainty. 

Combining these two results we think that a reahstic value of the total elec- 
tromagnetic correction is 



The error given in Eq.(31) is a generous one; it may be smaller. There is some 
uncertainty due to the neglect of higher orders in the chiral expansion and we 
have allowed the possibility of a larger uncertainty from the constants li, 
4. Since the DIRAC experiment [1] aims to deliver a result for r that will give 
a value of with an error of 5 %, the uncertainty in Aoc(/^o) is much smaller 
than the expected experimental error. 

This means that, if the DIRAC experiment achieves its proposed accuracy, it 
will be able to give a result for aQ^(/io) with an accuracy comparable with that 
of the present estimates from ChPT. Recent work of Colangelo, Gasser and 
Leutwyler [18], in which two- loop ChPT is combined with the Roy equations 
and phenomenological data from pion-pion scattering, gives a° — a^=0. 265(4). 
This value is in units /i^^ and is based on a hadronic pion mass of /^c- Ftom 
the two loop values in Table 1 we obtain 



Aoc()Uo) 



0.055(5). 



(31) 



a°(/Xo) - a^iiio) = a^inc) - a^il^c) - 0.0143 fm 



(32) 
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If we convert the Colangelo et al. value to fm and assume that Eq.(32) is also 
true for these results, we then obtain 

«oc(/^o) = -^{a'M - a^i^o)) = -0.170(3) fm (33) 



Ref.[18] gives an error for aQcil^c) of only 1.8 %. However, it needs to be re- 
membered that the calculations were made within the framework of an isospin 
invariant effective hadronic theory in which the pions have the mass fic and 
that the phenomenological data used have not been corrected in any way for 
the presence of the electromagnetic interaction. The same remarks also apply 
to the determination of the low energy constants needed to make numerical 
predictions from ChPT. The error of 1.8 % could therefore be an underesti- 
mate. 

The value aQc(/io) = —0.170(3) fm, combined with the result in Eq.(31), leads 
via Eq.(ll) to the value r = 2.9 fs, which agrees with the result of Ref.[14]. We 
regard this value as an indication of the lifetime that the DIRAC experiment 
is measuring. The important result of Rcf. [17] and the present paper is that 
the electromagnetic correction in Eq.(31) is now determined with sufficient 
precision for a measurement of r with an accuracy of 10 % to give a value of 
the hadronic quantity aoc(A*o) with an accuracy of not much more than 5 %. We 
have also emphasised that true hadronic information needs to be obtained with 
a pion mass /iq. In particular, this mass needs to appear in the formal results of 
hadronic ChPT. Also, the phenomenological hadronic data used to fix the low 
energy constants and, at a later stage, to refine the hadronic scattering lengths, 
need to be corrected for the effect of the electromagnetic interaction. This is an 
important long-term programme and until it is implemented our knowledge 
of the hadronic scattering lengths will remain limited. On the other hand, 
since the electromagnetic correction given in Eq.(31) is now well determined, 
the DIRAC experiment will give a clean value of aQ^{no), whose error comes 
almost entirely from the error in the experimental result. 
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